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Abstract 

The quasielastic longitudinal electromagnetic response function R L is stud- 
ied within the context of a model that extends our previous investigations of 
pionic correlations and currents. Four mesons are now employed (tt, p, a and 
u, via the Bonn potential) and the many-body dynamics are extended to the 
full random-phase approximation built upon a Hartree-Fock basis. Wherever 
possible the Lorentz covariance of the problem is respected. The first three 
energy-weighted moments of the reduced response are computed, namely, the 
zeroth moment (Coulomb sum rule), the first moment (related to the position 
of the quasielastic peak) and the second moment (related to the peak width or 
variance). We discuss how with a modest downward adjustment of the Fermi 
momentum it is possible to obtain the expected zeroth and second moments; 
this implies that the nuclear system has a lower density than that required 
when using the free relativistic Fermi gas and accordingly that the interaction 
effects are weaker than one might initially find. 
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1 Introduction 



In past work [[!], Q we have addressed the problem of setting up both the longitudinal 
{Rl) and the transverse {Rt) responses that enter in the expression of the inclusive 
cross section for inelastic scattering of electrons from nuclei. The inelasticity domain 
of concern here is that of the quasielastic peak (QEP). Our studies were based on the 
relativistic Fermi gas (RFG) model and the guiding principles of our investigation 
were relativistic covariance and global gauge invariance (fulfillment of the continuity 
equation) — indeed much effort was expended in adhering to these two fundamental 
constraints. In particular, the choice of the RFG model was prompted not only by 
the assumption that the nuclear surface is likely of secondary importance in very 
inelastic nuclear phenomena, but also by the recognition that the two principles 
referred to here are easier to implement in a translationally invariant system than 
in a finite nucleus. 

Moreover, in our previous work we resorted to the Hartree-Fock (HF) and anti- 
symmetrized random phase approximation (RPA) many-body frameworks, adding 
their contributions to get the nuclear responses. These procedures are correct when 
the treatment is limited to first-order contributions in a perturbative expansion, an 
approach that has some merit when the emphasis is being placed on achieving a con- 
sistent treatment of the forces and currents (the meson-exchange currents (MEC) 
are in fact usually expressed through first-order perturbative diagrams) and when 
the interaction between nucleons is carried solely by the pion whose associated cor- 
relations are weak enough to make a first-order treatment quite reliable. Indeed, 
comparison with the infinite-order calculations showed this to be the case, at least 
for momentum transfers higher than roughly 300 MeV/c ||. 

In the present work we extend the previous investigations in two directions by 

i) including, beyond the it, additional mesons (namely the p, a and u) in the 
description of the nuclear force, and 

ii) calculating the nuclear electromagnetic responses in the fully antisymmetrized 
RPA scheme evaluated in a HF basis. 

Clearly, item ii) is implied by i), since the forces carried by the p, a and u are too 
strong to be dealt with in first order. On the other hand, the requirement of global 
gauge invariance is now much harder to implement in a consistent manner and we 
leave this issue to be dealt with in a future paper. 

Here we address the specific problem of finding how the forces between the nu- 
cleons shape the quasielastic response in the range of momenta between 300 MeV/c 
and 1 GeV/c, confining ourselves to the charge response Rl (the transverse one, Rt, 
is treated in a companion paper). To assess the impact of the forces on Rl we focus 
not only on the longitudinal response itself, but on three derived quantities as well, 
namely the sum rule, the mean energy and the variance of the response (we shall 



2 



later provide precise definitions for these quantities). Note that these observables 
relate only to the many-body content of Rl, which should therefore be carefully 
disentangled from nucleonic aspects of the problem as we shall see in detail. 

In performing the present investigation an underlying theme has been the role 
played by the pion. Apart from the self-evident importance of the pion in the 
nuclear dynamics, our interest has also been prompted by the findings of ref. ||, 
where it has been shown that, owing to its isovector nature, the pion acting alone 
induces a remarkable enhancement of the longitudinal vector weak neutral current 
response. This in turn entails a huge effect on a certain observable related to the 
asymmetry as measured in parity-violating inclusive electron scattering. Our wish 
to explore whether or not pionic effects survive when other powerful components of 
the nucleon-nucleon interaction are switched on is therefore understandable. 

Significantly, we have found that the role of the pion, far from being "sup- 
pressed", is actually "enhanced" by the shorter-range part of the nuclear force 
through an interference mechanism hidden in the complexities of the charge response 
of the correlated RFG. On the other hand, it is well-known that pionic effects in the 
ground-state energy of a spin-isospin saturated system like the RFG or a closed-shell 
nucleus are suppressed, because the key ingredient, namely the Hartree contribution, 
is eliminated by the spin-isospin traces. 

It thus appears that in nuclear matter the ground state and the QEP involve 
different aspects of the dynamics, namely that the action of the pion in the ground 
state is hampered, while in the QEP it is enhanced. This is a long-standing question 
that does not appear to have received unambiguous experimental resolution; how- 
ever, a recent reanalysis || of polarized (p, n) experiments at LAMPF shows that 
the data are compatible with a softening of the spin-longitudinal response induced 
by the pion, as predicted many years ago f5|. 

Since the developments and results presented in the body of this work are quite 
involved, in the following paragraphs we summarize some of the salient findings 
before entering into those more detailed discussions. In sects. |2] and | we build 
the HF mean field utilizing as input the tt, p, a and uj pieces of the well-known 
relativistic Bonn potential ||. We show that the HF field is strong and induces 
a very substantial hardening of Rl- This arises mainly from the tadpole diagram 
via the cooperative action of the a and the u, with however an additional sizable 
contribution having the same sign and stemming from the p meson through the 
"oyster" (Fock) diagram. Worth noticing is that the momentum dependence of the 
tadpole diagram, and hence the hardening of the charge response, is of relativistic 
origin. 

As we shall see, we are able to account quite accurately for the action of the 
HF field on the nuclear response simply by shifting the energy u transferred to 
the nucleus by an amount Auj and by changing the bare nucleon mass into 
a dressed one, m* N . It turns out that Auj and m* N depend upon the momentum 
transfer q and the Fermi momentum hp. This result is the key for obtaining RPA 
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results in a HF basis without having to resort to quite heavy computations. The 
strength of the mean field up to the largest momentum considered here, q — 1 
GeV/c, is also reflected by the variance of the HF longitudinal response (essentially 
the half-maximum width) which is significantly larger than that occurring in the 
RFG model, but for the lowest momenta. 

In sect. [| we investigate the roles played by the ring and exchange correlations 
in Ri. We deal exactly with the rings, whereas the exchange diagrams are treated 
in the framework of the continuous fraction method truncated at the first iteration. 
This scheme actually provides exact results for zero-range forces and also turns out 
to be quite accurate for finite-range interactions as long as the momentum transfer 
is not too small (say q > 300 MeV/c) 

In the first part of sect. f|we consider only diagrams going forward in time (anti- 
symmetrized Tamm-Dancoff (TD) approximation) and, to assess the impact of the 
mean field on the TD results better, these are presented with fermions propagat- 
ing both freely and in the presence of a HF field. In the TD scheme (and as well 
in RPA) a struggle among the different contributions to the charge response takes 
place. Specifically the battle arises from the interplay between the direct (ring) and 
exchange diagrams and among the various mesons carrying the force. To unlock 
this complex situation it might help to recall that 

a) the momentum behaviour of the ring and exchange diagrams is rather similar, 
both of them decreasing with q and becoming much reduced for q > 500 
MeV/c; 

b) generally, the light (heavy) mesons are more effective in transmitting the in- 
teraction at small (large) q, providing of course that the coupling constants 
are large enough. 

Thus, considering first the ring diagrams, it is found that they indeed contribute 
sizably up to quite large transferred momenta and are obviously supported only by 
the a and u mesons. These dramatically reshape the free charge response in the 
isoscalar channel when acting alone, leading to the occurrence of collective modes 
lying above (the u) and below (the a) the energy range of the QEP that carry a 
large fraction of the charge strength (hence leading to a huge depletion of Rl inside 
the QEP). However, when acting together, the a and the uj interfere so strongly that 
the "ring" collective modes disappear and the only remnant of the <j-lo interplay is a 
still appreciable damping (but for the lowest energies) and flattening at 300 MeV/c 
(but not for larger momenta) of Rl. In addition a softening of Rl takes place - 
this is large at 300 MeV/c, still significant at 500 MeV/c and then fades away at 
higher momentum transfer in accord with a). On the other hand, in accord with b) 
the a wins over the u, as it also does for the ground-state binding energy. 

Turning to the exchange, it is contributed to by all of the mesons. Furthermore it 
opposes the action of the rings, thus leading in accord with b) to a hardening of Rl, 
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which shows up as a dramatic peak in the response at q = 300 MeV/c, with some 
shadow of it remaining even at 500 MeV/c. This peak corresponds to an "exchange" 
collective mode embedded in the particle-hole (ph) continuum that arises from the 
repeated exchange of the attractive n, p and a mesons, which are more powerful 
than the u at these momenta inside a ph ring. Of considerable significance in this 
connection is the role of the pion and, to a much lesser extent, of the rho: indeed, 
without the pion the exchange peak would not appear at all! 

Actually, when the ring and the exchange act together, the repulsive peak disap- 
pears, although some hardening of Rl remains at 300 and 500 MeV/c, thus showing 
the dominance of the exchange correlations over those from the ring diagrams. At 
larger momenta the ring and exchange diagrams, beyond being quite small, almost 
completely cancel each other. 

Summarizing: the TD collective effects in the charge channel arise basically 
from the exchange diagrams. They show up essentially through a hardening of the 
charge response at intermediate momenta (due to the attractive a, but particularly 
to the 7r), which persists up to about 500 MeV/c. However, as we shall see, the 
TD correlations appear to be somewhat overcome by the HF field: thus the TD 
longitudinal response evaluated in a HF basis, while being appreciably depleted 
at intermediate momenta, is actually shifted at higher frequencies for all of the 
transferred momenta considered by an amount which is rather constant but for the 
lowest momenta. 

In the second part of sect. f| the diagrams going backward in time are explored 
(antisymmetrized RPA). The ground-state correlations (gsc) thus introduced appear 
to strengthen the action of the attractive ring diagrams somewhat, while at the same 
time reducing the contribution of the repulsive exchange. Indeed the "exchange 
peak" at q — 300 MeV/c, referred to above in the TD scheme, now no longer 
shows up and is instead replaced by a "ring peak" , associated with the a meson and 
appearing in the direct channel. Yet the non-linearity of the RPA is such that, even 
in the presence of gsc, the exchange term still prevails over the direct one when they 
are combined together and what remains is actually an appreciable hardening. In 
addition to this, the RPA correlations at momenta up to about 500 MeV/c (and 
likewise the TD ones, but more markedly so) have a tendency to quench and flatten 
Rl, particularly when evaluated in a HF basis, namely they produce trends that 
appear to be supported by some experimental evidence 0-[[Ll|. At larger momenta, 
the contribution of the gsc essentially vanishes, as can also be inferred from the 
observation that here the sum rule is well fulfilled. 

The important role that the pion also plays in the RPA framework should be 
emphasized. Indeed, without it the ring peak would survive the confrontation with 
the exchange contribution: accordingly, the RPA response in a free basis, and in a 
HF one as well, would still display a peak on the low-energy side of the response. 
This would be entirely out of touch with reality, if only the a and the ui were to act, 
and still there in a scheme encompassing the p as well. Thus, the statement that 



5 



the RPA charge response requires the pion appears to be justified. 

In summary, the RPA correlations are found to yield an Rp that is appreciably 
less hardened than the one obtained in the TD scheme, but somewhat more quenched 
and with a slightly larger variance. This finding is in accord with the general result 
that the gsc lower (enhance) the response when the latter is hardened (softened) in 
order to comply with the energy- weighted sum rule. Finally, as in TD, the HF mean 
field still tends to play a prominent role within the RPA and again the role of the 
pion is crucial. 

From our analysis of the longitudinal response based on the HF TD and HF- 
RPA schemes with the Bonn potential as an input, we thus conclude that the corre- 
lations' effects are very substantial and to a large extent appear to disrupt the RFG 
predictions. One is accordingly led to address the question whether or not part of 
the interaction might be accounted for by "renormalizing" the only parameter char- 
acterizing our many-body system, namely the Fermi momentum fcp. Our choice, 
namely kp = 225 MeV/c, has been dictated by the analysis of the 12 C data where 
such a value turns out approximately to account over a range of momenta for the 
width of the longitudinal response when the latter is analyzed in terms of the RFG 
model. Since the width at half-maximum together with the area and position of 
the QEP are the central features of the inclusive charge response (or, in general, of 
any response function), we have searched within the scheme of our meson-correlated 
RFG for a kp capable of accounting for the experimental half- width or, which is 
roughly the same thing, of yielding the same half- width as the RFG with kp = 225 
MeV/c. In other words we have applied what can be defined as a principle of the 
invariance of the variance, the latter being of course related to the half-width of the 
response. The integrated response and peak position are then computed using this 
renormalized model. 

The results are discussed in sect. [5]. We have been successful in adhering to this 
principle, providing we let kp vary slightly with q. The values of kp necessary for 
achieving this turn out always to be lower than 225 MeV/c by about 10% and to 
be gently increasing with q (i. e., the larger the value of q, the smaller the required 
renormalization of the Fermi momentum). Because of the non-linear dependence on 
kp of the interaction effects in our model, the impact of the correlations on Rp is 
overall found to be appreciably diminished. Note that, due to the different ranges of 
the forces carried by the different quanta, the heavier mesons are those most affected 
by the procedure outlined above, whereas the pion emerges almost unscathed from 
the renormalization of kp. 

Whether the results obtained in the invariance of the variance scheme lead to 
closer contact with data or not is a question that remains to be answered in detail. 
This we have not done in the present work partly because we have not yet com- 
pleted the evaluation of the MEC contributions, which should of course be added 
before attempting in-depth comparisons with the data, but additionally because two 
particle-two hole (2p — 2h) excitations cannot be ignored in the QEP domain and 
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should be included as well in the analysis. 



2 The Hartree-Fock field 

Two diagrams contribute to the HF field, namely the direct (sometime referred to 
as tadpole) and the exchange (sometime referred to as oyster) diagrams displayed in 
Fig. H We have calculated their contribution to the nucleon self-energy utilizing the 
Bonn potential, which is derived through a non-relativistic expansion of the diagrams 
corresponding to the exchanges of the ir, p, a and u mesons (actually the a meson 
should be viewed as a fictitious entity corresponding to the exchange of a correlated 
pair of S'-wave pions). Here we have neglected the small r\ and 5 contributions in 
ref. || and for clarity retained only the largest pieces of the potential. It reads 



V = V* + V p + V a + V", 



where 



V 7 ' 



V 



V u 



{V$<T 1 - t T 2 +V]!S 1 2)T 1 -T2 
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mr. 



r 
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In the above S\ 2 is the usual tensor operator and 

A 2 - m 2 



r«(g) 



A 2 +g 2 



(6) 



the standard meson-nucleon form factor. Moreover, q = |q| is the momentum ex- 
changed between the two nucleons and P = |P| is the average of their relative in- 
coming and outgoing momenta (as illustrated in Fig. |6], P = (ki — h 2 + k' x — k' 2 )/4). 
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meson 


T 


.r 


m (MeV/c) 


A (MeV) 


S747T 


/74tt 


7T 


1 


0" 


138.03 


1300 


14.9 




P 


1 


i- 


769 


1300 


0.95 


0.95 x (6.1) 2 


O 





0+ 


550 


2000 


7.7823 




UJ 





1- 


782.6 


1500 


20 





Table 1: The quantum numbers, the masses, the cut-off parameters and the coupling 
constants of the mesons included in the Bonn potential. Note the two coupling 
constants characterizing the p. 

The P-dependent terms in the p, a and u potentials clearly correspond to a non- 
local interaction in coordinate space and their important impact on the dynam- 
ics of the system will be discussed later. We have omitted the spin-orbit terms 
(oc ((Ti + cr 2 ) • q x P) because they do not contribute to the responses in an infinite 
system. 

For the convenience of the reader we report in Table [I] the quantum numbers, 
the values of masses and cut-offs and the strengths of the couplings for all of the 
mesons entering into the potential. 

Note that in our approach, as in refs. [I], 0, the short-range physics is again 
embedded in phenomenological vertex form factors which cut off the nucleon-nucleon 
interaction in a spatial region of size ~ 1/A around each nucleon. Whether this 
is adequate or whether ladder diagram contributions to the responses should be 
explicitly accounted for in the QEP domain is unclear. This will need to be addressed 
in future work. 

In connection with the HF field, we recall that only the a and u contribute to 
the tadpole term. Introducing dimensionless parameters 

k F A a m a 

VF = , K = , fJ>a = , (7) 

m N m N m N 

and measuring both the momentum and the self-energy in units of the nucleon 
mass ((3 = k/rriN, = mjvS(/3)), the corresponding real and energy-independent 
self-energies read 

On the other hand all of the mesons contribute to the exchange (oyster) term and 
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the associated self-energies, which are again real and energy- independent, read 



where 



and 



+^[SM + SA-m (12) 



° ' /: ' [&(/?) + &(-/*)], (13) 



16 4vr 2 



(14) 



+^gfaM - « + 2A^ - A*)ln[(^ + (3f + A 2 



+ 2! (/?2 " & ~ ^ )ln[(?7F + ^ + ^ (15) 



We thus see from the above that the Hartree field stemming from the action of 
the cr and of the u mesons displays a term that is quadratic in the momentum. The 
origin of this lies in the pieces of Vq and Vq that are proportional to (q/2mN) 2 and 
(P/2m7v) 2 , namely, to the parameters characterizing the non-relativistic expansion 
of the Bonn potential. In this connection note that the average between the initial 
and final relative momenta of the interacting nucleons in the Hartree diagram is 

P = ^(k-k'), (16) 

k being the momentum of the propagating nucleon (particle or hole) and k' the 
momentum flowing in the bubble of Fig. 0. Accordingly, in principle P is unbound. 
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However, in practice, in order that the truncation of the non-relativistic expansion 
of the Bonn potential makes sense, both |P| and |q| should not exceed, roughly, the 
nucleon mass m^. 

Because of its quadratic dependence on the momentum the Hartree field can be 
accounted for by an effective mass (of course, but for the constant terms which enter 
into (|) and ©): 



m 



m N 



N 



1 + 



4:171% 



/4 




(17) 



where 



P 



2k F 



m 



2rf F 



3tt 2 JV 3vr 2 v ' 

is the uniform density of the Fermi system. The effective mass is displayed as a 
function of hp in Fig. |6], where the individual contributions to m* N arising from the 
a and the u are shown as well. It is worth noting that the a and the u cooperate 
in dressing the bare mass of the nucleon whereas, as is well-known, they fight each 
other in accounting for the binding energy of nuclear matter. Clearly, this outcome 
reflects the opposite sign of the constant terms and the same sign of the quadratic 
terms of the self-energies in (§) and (§). 

Turning to the Fock contribution, the nearly constant behaviour as a function of 
the momentum of the combined Fock self-energy of the a and u (see Fig. ||) is quite 
striking — note that a constant here produces no effect on the HF response of the 
system. On the other hand the Fock self-energy of the p is sizable: actually, if acting 
alone, being repulsive, it would harden the charge response with respect to the free 
case by about 20 MeV at q = 500 MeV/c for k F = 225 MeV/c. It clearly cannot be 
incorporated in a constant effective mass over a wide range of momenta because, far 
from being quadratic, the Fock self-energy of the p displays a rather cumbersome 
momentum dependence (see eq. [IT]). Finally, we recall that in the oyster diagram 
P = 0. 

In order to appreciate the dynamical impact of the HF field, let us now calculate 
the longitudinal electromagnetic response in the HF approximation. For this purpose 

H 



we start by recalling its well-known expression 



namely 



3n 2 A 



lc 3 



x5 



u — 



Q 



dk 

- qk 



2m N m N 



6(k F -k)8(\k + (i\-k F ) 

-S HF (|k + q|) + S HF (fc) 



(19) 



A being the nuclear mass number (here and in the following we take N = Z = A/2), 
Q the four momentum (u, q), fi{Q 2 ) the longitudinal 7NN vertex and E HF the HF 
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self-energy. It is interesting that (|19f) can be expressed as an integral over the 

U l 



effective mass according to [T3| 



with k = q/2m>N, £a = 3ir 2 A and 



7i z r]pKm N 



(20) 



9{7) F -2k) / dpp 



If 



+9(r] F - k) 6{2k - Vf ) r d(3(3\m* N {JAK 2 + {3 2 + AKy>) 

J2K-ri F J v 

+ {6{K-rj F ) I dp (3 + 6(ri F - k) 9(2k - r] F ) J dp p\m* N (yJ An 2 + P 2 + Any"). 

(21) 



In the above 



m* N (P) 



m N 



1 + 



1 dZ HF (P) 



P dp 



(22) 



is the momentum dependent effective mass of the nucleon, whereas y' and y" are 
roots of the equation 



1 



HF 



0. 



satisfying the constraints 

{rf F - Ak 2 - P 2 )/Ak <y'<P 



and 



P < y" < P- 



(23) 



(24) 



It is of importance to realize, in connection with the basic result fl2~T|), that in i?^ F 
(and also in the transverse HF response) only the particle (and not the hole) effective 
mass explicitly appears. 

Notice that the relativistic content of our approach is contained first of all in the 
relativistic scaling variable of the Fermi gas ip. For the latter we use, rather that 



the exact one flljj, the approximate but quite accurate relativistic expression 



Vf 



A(A+ 1) 



(25) 



with A = uj/2m^. Next, relativity is also embedded in the longitudinal 7NN vertex. 
For this we again employ the approximate but almost exact p] expression 



fl(Q 2 



;i + 2A) 



1 



1 + r 



[G^(r)] 2 + r[G^(r)] 2 ^(l-^ 2 ) 



(26) 



11 



q (MeV/c) 


A (10 2 ) 


B 


A (10 2 ) 


B 


m* N (p)/m N 


300 


-1.671 


0.246 


-0.961 


0.146 


0.77 


500 


-1.671 


0.246 


-0.424 


0.116 


0.81 


800 


-1.671 


0.246 


1.034 


0.085 


0.85 


1000 


-1.671 


0.246 


1.466 


0.081 


0.86 



Table 2: The dimensionless parameters characterizing the biparabolic nucleoli 
self-energy discussed in the text (see ( p7f) and (p8| )). Note the variation of the 
particle self-energy with the momentum transfer q, reported in the last column. 
The effective mass of the particles below the Fermi surface (holes) turns out to be 
m* N (h) = 0.67 m N . 



Ge* m being the proton and neutron (with isospin projections m t ) electric and mag- 
netic Sachs form factors and r = \Q 2 \/4m 2 N . 

Before displaying the results for R^ F , we now briefly outline an approach that 
greatly simplifies the RPA calculation of the nuclear responses based on a HF (rather 
than on a free) basis. This approach stems from the observation that the total HF 
nucleon self-energy induced by the n, p, a and u mesons, while not simply parabolic 
in the momentum, nevertheless lends itself to an interpolation with two parabo- 
las, one below (hole states) and one above (particle states) the Fermi momentum. 
Accordingly, we set 

S HF (/5) = A + B{3 2 } for p < rj F (27) 



and 



£ ((3) = A + B(3\ for (3 > rj F , (28) 



B and B being both positive. However, while A and B are truly constant, A and B 
are in fact momentum transfer dependent, since they are fixed by fitting the particle 
self-energy over the range of momenta required for the calculation of a given nuclear 
response. This range, of course, varies with q. As already noted, the procedure 
outlined here allows a quite accurate fit of the self-energy entering into the charge 
response. 

In Table |2| we give the values of the parameters appearing in (|27|) and (p8|) for 
q = 300, 500, 800 and 1000 MeV/c. 

Now with the self-energy ( |27H28|) the solution of eq. (|23|) reads 

V = kT+Tb ~ K + T+~2B + < 29 > 
with a = (A — A)/2 and b = (B — B)/2, and the constraint (|24]) (non-Pauli-blocked 
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region) is easily shown to require 




l + 2B)n 



Then, by expanding the square root in leading order and noticing that 
self-energy ( p8|) yields the expression 

m N 



(30) 



with the 



m N = (31) 

for the particle effective mass, we obtain for the lower limit of the momentum inte- 
gration over the effective mass in ( ^l|) the value 



0o 



1 



1 + 25 

in the non-Pauli-blocked region. In ([E 



(X* + a*)- 

K 



K 



(32) 



|) the star reminds us that the effective mass 
(pl|) should be used in the definition of the dimensionless energy and momentum 
transfer variables and a* = (1 + 2B)a. 

Finally, by inserting (^) in the expressions (|20| ) and ([21]) for the longitudinal 

response, it is immediately apparent that one gets for the HF response the same 
expression as for the free response but for the replacement of the bare nucleon mass 
m N with the dressed one m* N and for the shift 



A ^ A* + a* 



(33) 



of the energy. 

Identical results can be shown to hold in the Pauli-blocked region, although in 
this case, owing to the energy shift (|33|), the response obtained with the approxima- 
tion outlined here no longer exactly vanishes as A goes to zero as it should. On the 
other hand we should not trust the physics of the RFG in the Pauli-blocked region 
anyway. 

The above very useful and appealing results are crucially dependent upon whether 
the leading-order expansion of the square root in (|30| ) is sufficient or not in providing 
a good representation for the lower limit of the momentum integral entering into 
the expression of Rl- That this is indeed the case will be shown in the next section. 



3 The charge response in HF 

We have calculated Rl in the HF scheme with the self-energies discussed in sect. 
for q = 300, 500, 800 and 1000 MeV/c.The results are displayed in Fig. ^| where, in 
addition to the charge response of the free Fermi gas, three cases are shown as well, 
corresponding to the HF evaluation of Rl with 
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i) the true self-energy S HF ; 

ii) the self-energy fitted with a double parabola, as explained in the previous 



iii) the analytic method, which we recall amounts to dressing the nucleon mass 
and shifting the energy. 

Two features are immediately apparent from the figure. The first relates to the 
outstanding agreement between the results for Rl obtained with the three prescrip- 
tions referred to above. The expansion of the square root in ( |3~0D indeed appears to 
be warranted over a large span of momentum transfers. The second feature relates 
to the systematic hardening of the response induced by the HF mean field when 
compared with the free response. This appears to be substantial, is persistent as 
q grows and, as already pointed out, stems mostly from the quadratic piece of the 
Hartree self-energy of the a and the u. 

In this connection, it is worth recalling that for local static potentials the Hartree 
self-energy does not contribute to the response. It does so in our case because the 
potentials mediated by the a and the ui also depend on the average between the 
initial and final relative momenta of the two interacting nucleons, P, as seen from 
the non-relativistic reduction of the Bonn potential at order 0(l/mjv)- 

To gain further insight into the impact of the HF field on the nuclear responses 
and to pave the way for a similar analysis of the additional correlations we shall 
consider later, we now investigate the sum rule H, the mean energy A and the 
variance a of Rl- These quantities should relate to the many-body content only of 
the nuclear responses, which should accordingly be disentangled beforehand from 
the complete charge response. To achieve this we follow the pattern established in 
ref . || , where it has been shown within the context of the RFG that a reduced charge 
response S l (k,X) can be introduced, having an energy integral (the Coulomb sum 
rule) equal to one for k > rjp- It is defined as follows: 



section; 



S l (k,\) 



vlRl 



(34) 



where 




(35) 



is the usual leptonic factor and 




(36) 



( KVF /2£ F )(d^/d\) ■ 



In the above, £p 




dip 
dA 



K ^i + ew 2 /2 ri + 2A + o^ 2 



(37) 



r [ 1 + A + CfiP 2 
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K ' 1+2A Wi 



and 



with 14 



W 2 
A 



1 + r 



7] F T \ 1 + A 



X L = v L —(G 2 E + W 2 A) 

T 



[G 2 E + rG 2 M ], 

(4 + e F T + T 2 ) + X(e F + T) + A s 



:i + r) 



and 



T = max ^ (e^ — 2A) , kJ 1 H A f . 



(38) 
(39) 

(40) 
(41) 

(42) 



As in (^) M reminds us that in taking the ratio (0) a term with M = Z should 
be added to a term with M = N in both the numerator and denominator. 

In terms of the reduced response (@), the Coulomb sum rule, the mean energy 
and the variance, namely related to the zeroth, first and second moments of the 
charge response, then read as follows: 



duS L {K, A), 
du A S l (k, A), 



and 



a 



du S l (k, A) [A — X L ] 2 



(43) 
(44) 

(45) 



In discussing our results, we recall that the HF sum rule should be identical to the 
RFG one. This is indeed seen to occur in Fig. ^a, where the momentum behaviour of 
and Spp G is displayed. In Figs. ||b and |^c we display the momentum behaviour 
of Arf, Arfg and ctrf, <5rfg, respectively, for the charge response. 

It clearly appears from Fig. []b that up to the largest momenta considered the 
HF field hardens the charge response as previously mentioned. With respect to a, 
it can be shown for k > rj F in the RFG framework that it is related to the width 
Au 



m7v(4K^/v / l + 4k 2 ) of the response region according to |0| 

1 qk F 



a 



2mjv\/5 

1 HT] F 



mjf + q 2 



+ 0{k 



VI VI + 4k 2 

—= + 0{rf F 

4V5 m N 



+ 0( Vf ) 



(46) 
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Remarkably, in Fig. |6]c one observes that in spite of the dramatic hardening of the 
charge response induced by the HF field, which grows with q, the difference between 
the variance of P^ F and of P^ FG is quite small at low momenta (~ 5%), then it 
increases with momentum transfer up to about q = 500 MeV/c (~ 20%), to decrease 
again for larger momenta where the difference stabilizes around a value of 10%. 

Thus, notwithstanding the violation of the energy-weighted and energy- squared 
weighted sum rules induced by the HF mean field, which is in accord with a general 
theorem on the contributions of the Feynman diagrams to the sum rules yet 



it appears that these violations tend to compensate each other in arriving at ctrf- 
It is as if the system reacts against those interactions that attempt to change the 
width of the charge response or, equivalently, the Fermi momentum kp, which in 
turn fixes the density for a homogeneous system. 



4 The RPA longitudinal response 

In this section we explore the influence on the charge response of correlations between 
nucleons acting beyond the mean field. It might be useful in this connection to 



relate the ingredients used in the present work to Feshbach's scheme |jT7 |, which 
splits the Hilbert space of the nuclear Hamiltonian into a P-sector, where the slow 
nuclear motion occurs, and a Q-sector, where instead the fast motion associated 
with short-range correlations between nucleons takes place. In this framework the 
HF field and the RPA correlations are naturally included in the P-sector, whereas the 
2p — 2h, 3p — 3h, etc., types of excitations mainly belong to the Q-sector. Indeed the 
occurrence of the latter requires a large amount of energy which is mostly supplied by 
the violent nucleon-nucleon interaction at short range and short-range encounters 
between nucleons are presumably comparatively rare for kinematics in the QEP 
region. Therefore, our aim will be to calculate the polarization propagator H(Q) in 
RPA on a HF basis while ignoring the Q-sector contributions in the present work. 

For this purpose, we start by recalling the link between H(Q), which for a homo- 
geneous system only depends upon the two variables |q| and cu, and the ph Green's 
function G ph , namely 

n(Q) = l Iwi~J Wj^ GPh{k + p + <?, p) ( 4? ) 

(spin-isospin indices are here neglected for simplicity). Therefore, in order to get n, 
one should first solve the Galitskii-Migdal equation [0 for G ph which reads in the 
RPA framework 

G ph (k + Q i k;p + Q,p) = -G(p + Q) G(p) (2n)%k - p) 
+iG(k + Q)G(k) J AL[u ph ^{Q)-U ph ^ c \k-t)] G ph (t + Q,t;p + Q,p), 



(4* 
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where the single-particle Green's functions G are meant to describe the fermion 
propagation in the HF field. 

Clearly, the zeroth-order iteration of fl48|) , when inserted into (|47D , just yields 
the polarization propagator for the non- interacting relativistic Fermi gas U°(Q). 
The first iteration of (f48|) , on the other hand, provides, via the direct ph matrix 
element of the interaction JJ ph ' An ', the first ring (Fig. |]a) of the ring expansion 
and, via the exchange ph matrix element JJ P ,e3a ^, the first (Fig. |6|b) of the infinite 
series of exchange terms. When the two expansions are grouped together, the fully 
antisymmetrized n RPA on a HF basis is obtained. Note that the two first-order 
terms of the RPA series already embody contributions going both forward (TD) 
and backward (gsc) in time. 

In this paper, rather then directly solving eq. (0), we resort to the equivalent 
continuous fraction representation of Il RPA (q, uj) (T^, [2D], |T7|. A detailed analysis of 
the method will be given in ref. ||21|| , where, in particular, it will be shown how the 



ground-state correlations can be incorporated in the treatment (see also ref. p2 
Here, to help in understanding the present calculation, we merely recall the basic 
formulas. 



In ref. [19 it had been shown that the polarization propagator in the TD ap- 



proximation and in the first-order continuous fraction expansion is 

n TD (q u ) ^ n fi(q^) (49) 



where the subscript R means that only the retarded component of the free propa- 
gator is included and where 

4 1} (q,u;) = -/ ^ / ^G R (A; 1 + Q)G^(fc 1 )G^ 2 + g)G^(A; 2 )\/(k 1 -k 2 ). 

(50) 

Moreover V(q) and V(ki — k 2 ) represent, respectively, the direct and the exchange 
ph matrix elements of the Bonn potential in the appropriate spin-isospin channel. 
These formulas can be generalized to encompass the ground-state correlations as 
well simply by inserting into them the full free propagator G° and n° (i. e., by 
including both retarded and advanced components). Accordingly the polarization 
propagator will read 

n RPA r n , a n °(q^) (ci \ 

i-\/(q)n (q^)-[na)(q,a;)/no(q,^)]' 101 } 

with 

n (1) (q^) = -j^ij ^G (fc 1 + Q)G°(fc 1 )G' (A; 2 +g)G (fc 2 )K(k 1 -k 2 ). (52) 
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The HF correlations should be included, in principle, using in eq. ([51]) I1 HF in place of 
n° and in eq. fl5"2]) G HF in place of G° (and analogously for the TD expressions). To 
simplify the calculations, we shall make use of the prescription discussed in sect. |2|, 
which amounts to replacing the bare nucleon mass with the dressed one and to 
shifting the energy (see the discussion above eq. Q33D). 



In the following we discuss the results that formulas ([491) and (|5T| ) yield with the 
Bonn potential. We split our analysis into two parts: in the first one (subsection 
[4.1|) we deal with the charge response in the TD approximation both in the free 
(i?™' frce ) and in the HF (_R™' HF ) basis; in the second one (subsection f4.2| ) we 
explore the charge response in the full RPA, again in the free (i?^ PA,free ) and in the 
HF (Rf A ' HF ) basis. 

4.1 Tamm-Dancoff correlations 

Note that in obtaining free and Srfo HF we have not used the normalizing factor 
(|36l), but rather the "relativized" version of it ||, namely 



x Lm=M^rizn»t> ( 53 ) 



ktjf dip/d\' 

since the TD charge responses in the free and HF bases are only approximately, 
although quite accurately, relativistic in the sense previously discussed in sect. |2|. In 
(|53|) ip is meant to be given by (^5)) and /| by (|26]) for the proton (Af = Z) and for 
the neutron (Af = N) , respectively. 

The reduced longitudinal responses S^ Diree and St DH F' like S# F , should yield 
the same sum rule as the non- interacting RFG (TB|. In Fig. ^a this is seen indeed 



to occur with a good accuracy. In connection with the mean energy A and the 
variance a, it is of interest to recognize in Figs. |6|b and c that the TD correlations 
alone (i?™' free ) do not modify either of these two observables with respect to the 
predictions of the non-interacting Fermi gas, although in principle they do not have 
to fulfill the energy-weighted and energy-squared weighted sum rules. This outcome 
is reasonably clear for the variance as seen in Figs. |6|-P, where Rl is displayed in HF, 
TD and HF-TD approximations at q = 300, 500, 800 and 1000 MeV/c for k F = 225 
MeV/c, whereas it is somewhat less evident that there is a near identity of the mean 
energy of the RFG and of the free basis TD responses at low momenta. 

To provide an orientation for understanding both the moments of the charge re- 
sponse and Rl itself over a wide range of momentum transfers we start by discussing 
the momentum dependence of the ring and exchange diagrams. In the former case 
such a dependence is mainly dictated by the direct ph matrix element of the Bonn 
potential, whose sole dependence is upon the momentum transfer q (and this is sub- 
stantial only in the isoscalar channel — in the isovector one the ph matrix element 
of the p is negligible): 

V(q) = V°{q) + V£{q) 
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2 



9 i 9 -7<T <7 \^ / 9 , 9 

g z + m^j q z + 

^ r '(g)-^-T+^(g)- 1 



q z + + m z 



(54) 



This formula most clearly shows the battle going on in the ring framework between 
the a and the u, and won by the former, as is apparent from Fig. ^| where V(q) is 
displayed. There it is seen that for momenta up to about 500 MeV/c a substantial 
attraction occurs that induces a softening of /^ in s> free _ This softening is evident in 
panel (b) of Figs. ^| and |6|. In contrast to this, for larger momenta V(q) is no longer 
strong enough to organize much collectivity in the nuclear motion and thus fails to 
affect Rl significantly (see panel (b) of Fig. |6| and ft should be remarked that 



strictly speaking formula (|54j) is not valid because the P-dependence of V has been 
ignored. On the other hand this seems to be warranted, since in the direct (ring) 
channel P = (k — k')/2, k and k' being the hole momenta of two subsequent rings. 
Thus P is quite small, being bound by kp. 

Turning to the exchange we first observe that in this channel P = q/2 and 
therefore the corresponding terms in the Bonn potential cannot now be neglected. 
Moreover, the interaction lines in the exchange diagram do not depend upon q. Nev- 
ertheless, for short-range interactions, the direct and the exchange matrix elements 
are actually rather similar: thus the effect of the exchange and likewise the ring dia- 
grams is no longer perceived at large q, the reason however stemming not only from 
the interplay among the different mesons, but also from the competition between 
the isoscalar (r = 0) and the isovector (r = 1) contributions, as will be illustrated 
in the following. 

In addition to the momentum dependence it is also of importance to recall the 
role played by the mass of the quanta exchanged by the nucleons. Indeed, the 
larger the mass, the shorter the range of the force and the slower the decrease with 
momentum of the interaction. 

On the basis of the above considerations the following comments are appropriate 
concerning the results displayed in the Figs. 

i) As already anticipated in the Introduction the forward-going ring diagrams, 
contributed to only by the a and u> mesons, dramatically affect the charge response, 
however in "opposite" directions. Indeed, they lead to the occurrence of isoscalar 
collective modes above (the u) and below (the a) the QEP energy range, thus 
strongly depleting the charge response itself in the QEP domain. Actually, the a 
taken alone would truly induce a phase transition in the RFC However, when the a 
and the to are considered together in the ring framework, their effect is substantial 
but not dramatic, yielding for q < 500 MeV/c an i?™' free that is quenched, flattened 

and softened with respect to i?^ FG . Then these effects fade away and at roughly 800 
MeV/c little is left of them (see the panel (b) in Figs. |§-|6|). This prevailing of the a 
over the u in the forward-going ring diagrams is in accord with the above discussion 
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and clearly goes in parallel with the analogous effect occurring in the binding energy 
of nuclear matter. 

ii) All of the mesons contribute to the exchange diagrams. Here the a and uj 
counteract the attraction and repulsion, respectively, that they induce through the 
rings, thus providing an exchange-repulsion in the case of the a and an exchange- 
attraction in the case of the uj (in the exchange an additional minus sign occurs). 
Furthermore, an important contribution that is substantial up to about 500 MeV/c, 
particularly in the isoscalar channel, is provided by the pion. 

Where the effects on the charge response are concerned, it should be observed 
that at q — 300 MeV/c an exchange repulsion prevails. It leads to the occurrence 
of a dramatic peak in ^ xch ' free ( see pjg_ already referred to in the Introduction. 
Fig. |6] is instrumental in grasping the origin of the peak. Specifically, in Fig. |]a the 
total exchange response ^£ Xch ' free j s displayed and compared with the one obtained 
using a force mediated only by the a and uj (or the a, uj and p). The pion's role 
appears to be truly stunning. Without it no peak would show up in the exchange 
response. Is this effect associated with the direct pionic contribution to i?™' free or 
with the interference that the pion is experiencing with the other mesons? 

The answer to this question can be read in Fig. ^|b where the i?™' free one would 
get with all of the contributions, except for the one arising from the interference of 
the pion with the other mesons, is compared with the full TD charge response in a 
free basis. The large difference between the two curves is then to be ascribed entirely 
to the existence in the exchange diagrams of the interferences between the pion and 
the other mesons. Note however that in the same figure the contribution AR^ due 
to the pion alone is also displayed and seen to be quite out of phase with that of 
the interference. Accordingly, the net pionic contribution to i?™' free is no longer as 
substantial as in the case of the pure exchange response. Yet it is sufficient, in spite 
of the presence of the large and attractive ring term, to bring about the hardening 
of the TD response seen in Fig. |6|d. 

For q = 500 MeV/c the peak is essentially gone, since the influence of the ir and 
the a on the one hand and of the uj on the other more or less balance each other, 
although some hardening of the response still takes place. It is worth observing that, 
in contrast to the nearly exact cancellation of their Fock contribution to the mean 
field, the a and uj do affect the charge response via the exchange term and that this 
arises mainly from the a, which acts in the same direction as the pion in shifting 
the maximum of Rl to larger energies for momenta up to about q = 500 MeV/c. 
For larger q the exchange becomes insignificant. This fading away is due not only to 
the compensation of the contributions of the different mesons (but for the p, which 
is very small by itself at large momenta), but also through an isoscalar-isovector 
compensation as illustrated in Fig. | for q = 1 GeV/c. In particular, there one sees 
that an appreciable pionic isoscalar exchange contribution still occurs in spite of the 
large momentum. 

An additional comment may be made about the roles of the it and p in the 
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exchange diagram. We first recall that the tensor force cannot act in the charge 
channel 0. Accordingly, the ir and p will be felt in only through the central 
part of the interaction they carry. For the 7r in particular, it is the short-range 
component (the (^-component for pointlike nucleons) which mainly comes into play, 
yielding, as we have already seen, a crucial repulsion at moderate momenta (indeed 
the pion is strongly attractive in the direct spin- longitudinal channel). In the case 
of the p, the finite-range component of the central interaction is stronger than for 
the pion: since it has an opposite sign to the zero-range component, a net, but much 
weaker, attraction mediated by the p thus results in the exchange channel. 

In summary for the exchange contributions, we find the following: the pattern 
of the total TD exchange contribution to Rl is quite complex — it turns out to be 
repulsive through the decisive action of the it that develops mainly in the isoscalar 
channel, which tilts the competition between the u and the a in favour of the latter 
for q less than about 500 MeV/c. Because of this repulsion an "exchange collective 
mode" appears as a marked peak on the high-energy side of the response for q = 300 
MeV/c. For larger q, the peak no longer appears and the exchange contribution 
becomes very small. 

iii) We now come back to the delicate balance taking place in the TD framework 
between the ring and exchange contributions. Indeed, when evaluated in a free 
basis, the trend of these is to cancel one another. The cancellation however becomes 
complete only around q = 800-1000 MeV/c. On the other hand for lower momenta, 
as we have seen (see panel (d) of Fig. |6|), a hardening of the charge response occurs, 
signaling the dominance of the exchange in the TD correlations that arises crucially 
from the 7r and a (since this is counteracting the u contribution). 

iv) We have previously mentioned the strength of the HF mean field arising 
from the Bonn potential. A question to be asked is then: what happens to the TD 
correlations when framed in the HF basis? 

Here, as before, the answer should distinguish between the range of momenta 
below and above 500 MeV/c. In the former instance (panel (d) of Figs. |6] and ||) 
a hardening, quenching and some flattening of Rl is clearly apparent. As in the 
case of the free basis it arises from the exchange diagrams via the 7r and a. For q 
above about 500 MeV/c an almost pure HF response shows up instead, which arises 
because of a substantial weakening of the force due to a compensation taking place 
among both the various mesons and TD diagrams. 

4.2 Ground-state correlations 

We treat the gsc in this separate section because of their subtle intertwining with 
the forward-going amplitudes, which significantly affects the interplay between ring 
and exchange discussed previously. 

From a glance at Figs. ^|-^| it is immediately apparent that the gsc strongly 
enhance the impact of the ring diagrams on the charge response, at the same time 
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damping the influence of the exchange contributions. More specifically, as already 
pointed out, the gsc go in the direction to suppress a hardened and enhance a 
softened R L and thus to respect the energy-weighted sum rule, as they are required 
to do. It is of relevance that the two effects, when they exist (namely for q < 500 

MeV/c), have a common origin. Indeed, we have seen in Fig. ^ that the direct 
ph interaction is always negative. On the other hand, neglecting for a moment 



the exchange term, it clearly follows from (51) that an enhancement of the charge 
response should be expected if 

Ren°(q,^)^^-. (55) 

In Fig. | the two sides of this equation are displayed versus u for q = 300 MeV/c. It 
is satisfying to realize that no solution of (|55|) occurs if only the retarded component 
of Rell is kept, whereas a solution exists when the advanced component (gsc) is 
included as well. Hence the peak on the low-energy side of Fig. |6|b and the absence 
of the peak in Fig. ||b are explained. Note that the existence or absence of a solution 
of eq. ( |55| ) does not depend upon the basis (either free or HF). Note also that the 
peak has a width, since it is embedded in the ph continuum. 

Furthermore, in the exchange channel the interaction turns from attractive to 
repulsive. The equation to be considered in this case reads 



Re 



n°(q,w) 



1. (56) 



Although eq. ([56]) can be solved numerically with and without the advanced com- 
ponents, it is of significance to realize that to the extent that 

n«(q,o;)«-y(q)[n (q,a;)] 2 (57) 

represents an acceptable approximation (which is the case for the present qualita- 
tive discussion), then instead of (|56|) one can consider (|55l), except that now the 
interaction is positive. A glance at Fig. || is then sufficient for concluding that gsc 
can prevent the occurrence of collective effects in the charge response for a repulsive 
interaction. Indeed, this is seen to be the case by comparing panels (c) of Figs. [] 
and H Moreover, by comparing Figs. |6| to |6| with the corresponding ones in the TD 
approximation, it is also apparent that the gsc are no longer really operative at large 
momenta. Thus, as in the TD case, in RPA we witness a vanishing of the impact of 
the correlations for increasing momenta and for the same reasons, namely, the fight 
between the a and u in the ring channel, the fight of the 7r and a against the u, 
plus the isoscalar-isovector cancellation in the exchange channel. 

The vanishing of the gsc at large q is clearly seen as well in panel (a) of Fig. ^ 
and indeed the Coulomb sum rule appears to be satisfied almost perfectly at q = 1 
GeV/c. At lower values of q one observes in the same figure that the gsc violate the 
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Coulomb sum rule, but more markedly so when they are evaluated in a HF basis 
rather than in the free one. On the other hand, from panels (b) and (c) of Fig. |6|, 
the mean energy and the variance of the response are seen to be essentially the same 
as those obtained in the TD framework. 

We should now comment on the critical role played by the pion in the RPA 
response. Notably, in RPA a situation arises that is symmetrical, but inverse to 
the one occurring in the TD. Indeed, the peak that shows up in Fig. ^b for the 
ring approximation (now containing components going both forward and backward 
in time) would still be present in both ^ PA ' free anc l _R^ PA ' HF ; except for the role 
played by the pion. This is clearly seen to happen in Fig. ^| where, as in Fig. |6|b 
but now in the RPA, the response one would get with all of the contributions added 
to the free one, except for the interference of the pion itself with all of the other 
mesons, is displayed and compared to the full response at q — 300 MeV/c. A broad 
peak at low energy sustained by the attraction provided by the a via the rings is 
still clearly present. In the same figure the ir-aup interference contribution is also 
shown together with the direct pionic contribution AR n . 

In contrast with the TD case, one sees that although these two terms are still 
out of phase, the dominance of the interference contribution is evident. One can 
thus conclude that at moderate momenta the pion is truly important in shaping the 
RPA charge response of the RFG and that it does so by interfering with the other 
mesons. It is accordingly true that the a and u, which are more strongly coupled 
to the nucleon than the pion, rather than washing out the impact of the latter on 
the charge response, are actually effective in enhancing it. 

5 Invariance of the variance 

The variance of the longitudinal response, which is related to the width at half- 
maximum of the response as shown by the RFG result (0), is defined in terms 
of the first and second moments of Rl- As already stressed it represents one of 
the basic attributes of the charge response and, in fact, any theory aiming to deal 
satisfactorily with the latter might be expected to account for this observable before 
attempting to reproduce finer details of Rl- 

Actually when experimental data are analyzed in terms of the RFG model, the 
Fermi momentum /jp is usually fixed in such a way as to reproduce the experimental 
width. Following this train of thought one could therefore search for a value of kp 
in the HF-RPA framework to give the same variance as that yielded by the RFG 
with kp = 225 MeV/c, the value that produces reasonable agreement for the width 
of Rl when compared with experiment for 12 C at q < 500 MeV/c. Note that fcp in 

the interacting model is expected to be different from the value obtained using the 
free RFG. We have performed such a search obtaining in the HF-RPA approach the 
values k F = 196, 201, 206 and 208 MeV/c at q = 300, 500, 800 and 1000 MeV/c. 
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In so doing we have arrived at a model of the longitudinal response that for large 
momenta produces the same zeroth moment (the sum rule, see Fig. |6]a) and, by 
construction, the same second moment or variance as the RFG at any q. 

Thus we see that the essential requirement of this procedure where the variance 
is maintained is to have a lower value for hp in the interacting problem than one has 
in the free RFG model. Correspondingly, one has a lower density in the former than 
in the latter. Over the range of momentum transfers explored — the range where 
the models may be expected to be most applicable — our results show that the most 
important feature is a renormalization downwards of kp by roughly 10% from the 
free RFG value. The residual g-dependence over the range 300 < q < 1000 MeV/c 
is quite mild. Indeed we have no reason to believe that it should be present at all, 
as its purpose here has been simply to enforce the principle of the invariance of the 
variance exactly. When in future work we attempt to make detailed comparisons 
with experiment (i.e., after the missing ingredients discussed in the Introduction 
have been incorporated in our model) the starting point for application of the present 
interacting model to 12 C will likely be with hp fixed to around 200 MeV/c. An 
immediate consequence of this statement is that the results presented here using 
kp = 225 MeV/c and the HF-RPA model will constitute our starting point when 
discussing medium-heavy nuclei, since the corresponding free RFG values of kp are 
in the range 267-252 MeV/c for q = 500-1000 MeV/c. 

Returning to Fig. |6|, we see that the basic global element of distinction between 
the HF-RPA and RFG models accordingly rests on the first moment of the response 
(i. e., the mean energy or energy- weighted sum rule). In Fig. ||b A L is indeed 
seen to be larger in HF-RPA than in the RFG almost uniformly in q (although, by 
comparing with Fig. |)]b, the gsc are seen to quench somewhat the difference between 
the correlated and free mean energies). 

To expand further on this issue we display in Fig. ||the g-behaviour of the shift e L 
of the maximum of the reduced charge response S l (k, A) with respect to that found 
for the RFG, \Q 2 \/2rriN. Indeed it appears that the HF-RPA shifts are somewhat 
lower than those of the HF-TD model, the difference disappearing at q ~ 1 GeV/c, 
in agreement with the findings of the previous section where the gsc contribution 
was shown to become insignificant somewhere in between 800 and 1000 MeV/c. In 
the same figure we also display the HF-RPA and the HF-TD predictions for e L for 
fixed kp = 225 MeV/c. The difference between the results in the two schemes clearly 
displays the effect related to the renormalization of kp — in fact the shift of the 
charge response is seen to be reduced by at least 35%, the larger reduction occurring 
at the lower momenta. 

The general trend seen here is an increase of e L with q up to around 800 MeV/c 
after which it begins to fall. The fact that e L > is borne out by most experimental 
measurements, although the trend with q is hard to be sure of, given the current 



state of the data — in fact little shift is seen at all in some cases 11 , while others 



show a g-behaviour that is quite compatible with the results we have obtained 
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at this (early) stage in our model-building. One thing is clear, however: the size 
of the shift seen in Fig. |6] is larger than that seen in the experimental studies to 
date. Clearly the ingredients that are still missing from our analysis should be 
incorporated before serious comparisons with experiemnt are made and clearly also, 
given the rather confused situation that still obtains on the experimental side, new 
measurements are essential. 

The HF-RPA charge response is shown in Fig. ^ over the usual range of momenta 
both for hp fixed at 225 MeV/c (corresponding to medium-heavy nuclei) and for 
hp varying with q according to the principle of the invariance of the variance (and 
roughly representing 12 C, see above). The overall quenching of the interaction's 
effects on the longitudinal charge response is clear, as is its evolution with q. Note 
that the quenching of the contribution to the longitudinal response of the individual 
mesons is g-dependent as well, an outcome possibly related to the different ranges of 
the interaction carried by the different quanta. As an example to illustrate this point, 
in Fig. ^| Rl is displayed for q = 300 MeV/c along with the individual contributions 
of the 7r and o . Here the renormalization of hp is the largest and accordingly the 
suppression of the interaction the strongest. One sees that the pionic contribution 
to Rl is significantly less affected by the change in hp than the one arising from the 
a. 

6 Conclusions 

This paper represents a further investigation (among many) of the nuclear charge 
response function. It therefore seems appropriate to comment briefly on its moti- 
vations. In addition, while most of our basic conclusions were anticipated in the 
Introduction, there are nevertheless a few that bear amplification. 

Clearly one of our motivations in the present work has been to explore Rl over 
a wider span of momenta than is usually done, having in mind possible future 
experiments to be carried out for example at M. I. T. /Bates, Mainz and CEBAF. In 
order to accomplish this task fulfillment of Lorentz covariance is of major importance 
and indeed much attention has been paid by us in this connection. Furthermore, 
in order to provide reliable predictions for Rl at momenta up to as large as 1 
GeV/c, a reliable force should be utilized. We have chosen the Bonn potential 
to be representative of modern approaches in this range of energies. Of course 
in calculating the longitudinal response the harmony between currents and forces 
must be preserved. Thus the predictions on R L offered in this paper should first 
of all be supplemented with MEC contributions before detailed comparisons with 
the available experimental data are attempted. This is a task we are currently 
undertaking. 

We have seen that the application of the principle of the invariance of the variance 
has led us to calculate the longitudinal response at values of k F somewhat lower than 
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those usually considered, the latter having been determined within the context of 
the free RFG model which should be expected to yield different results than our 
HF-RPA model. The consequences of this renormalization to lower values of kp 
and hence to lower densities for those contributions to Rl that we have neglected 
in the present analysis, namely the 2p — 2h excitations, the MEC effects and effects 
from ladder diagrams, remain interesting issues to be explored further. All of these 
scale more rapidly by extra powers of kp than the HF, TD and RPA contributions 
considered by us in the present paper. Accordingly, if we stick to the requirement of 
reproducing the variance of the experimental response within the HF-RPA scheme, 
their contributions to Rl are likely to be substantially suppressed. With respect 
to this principle it should also be observed that were our HF-RPA theory to be 
able to account for the mean energy (i. e., the first moment of the experimental 
response) then it would ipso facto account as well for the second moment. This could 
already bring the HF-RPA framework into rather close touch with experimental data 
and, as a consequence, less room would be left for the contribution of the MEC, 
2p — 2h and ladder diagrams. Our starting point in future work aimed at exploring 
these contributions will be to take kp = 200 MeV/c for 12 C and 225 MeV/c for 
medium-heavy nuclei. It remains to be seen what the ultimate predictions will be 
for the moments of the reduced charge response (including moments higher than 
two that characterize the shape of the response) and how these do or do not agree 
with experiment. Such detailed comparisons will be attempted once the missing 
ingredients have been incorporated in our extended model. 

Moreover, clearly even at this stage in our analysis we see that it is very desirable 
to have the experimental situation concerning the longitudinal response clarified, 
especially given its present status where it is characterized by conflicting data, since 
only then can definitive tests be made of our predictions. 

Be it as it may, our analysis shows that, even in the limited HF-TD and HF-RPA 
schemes, the major reshaping of the charge response with respect to the predictions 
of the pure RFG model arises from a rather complex intertwining among different 
forward- and backward-going diagrams and among the different mesons carrying the 
NN interaction. In our view it is of significance that there appears to be good reason 
to believe that, of the latter, the pion ultimately plays the crucial role in yielding 
the shape and the norm of the charge response of the correlated RFG. As illustrated 
in the present paper this is accomplished through quite subtle aspects of the nuclear 
many-body problem. For example, it is clear that the interference occurring between 
the pion and the other mesons acquires a quite different magnitude according to 
whether it occurs only in forward-going processes (TD) or in the ones going backward 
in time as well (RPA). It appears that precisely these interference terms significantly 
enhance the role played by the pion in the charge response. This also suggests that 
the effect of the enhancement of the weak neutral current longitudinal response that 
we predicted || in the framework of a purely pionic model will survive even in the 
presence of other mesons. 
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A relevant question to be addressed for the longitudinal response obtained in the 
present paper within the context of the RPA-HF model concerns its sensitivity to the 
parameters that characterize the pieces of the Bonn potential that we have employed. 
There are nine such parameters altogether, namely five coupling constants (two for 
the p) and four cutoffs, i.e., when the masses are taken as fixed which might be 
questionable in the case of the a-meson. 

The stability of the mean field with respect to moderate variations of the param- 
eters appears to be well-established: we have indeed seen that the mean field arises 
mainly from the quadratic terms in the self-energy of the a and uj and that these 
have the "same sign". 

The situation concerning the RPA correlations is different in that there a clearcut 
statement cannot be made. In the direct channel the competition between the a and 
uj mesons, both in forward- and backward-going rings, leads to a cancellation of the 
associated contributions which can be expected to be affected by variations of the 
force parameters. This, of course, is an effect analogous to the one that occurs for 
the ground-state energy. In the exchange channel the o-uj competition is altered by 
the action of the n and p, especially the former. Therefore in this case the sensitivity 
to variation of the force parameters will stem not only from contributions that tend 
to have delicate cancellations, but also from the non-linearity in kp contained in the 
pionic interferences with the a and uj. 

We are presently carrying out the task of assessing the extent of these sensi- 
tivities to details of the input potential model, despite the complexity of such an 
undertaking. 
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Figure captions 



Fig. 1 The two Feynman diagrams contributing in first order to the self-energy: (a) 
tadpole (direct), (b) oyster (exchange). 

Fig. 2 Definition of the momenta entering the Bonn potential: accordingly, one has 
q = ki - k; = k' 2 - k 2 and P = (k x - k 2 + k; - k' 2 ) /4. 

Fig. 3 The effective mass stemming from the Hartree self-energy of the a (dashed 
line) and of the uj (dot-dashed line) versus the Fermi momentum kp. The 
solid line represents the Hartree effective mass arising from the total a + uj 
exchange. 

Fig. 4 The Fock self-energy of the a (dashed line) and of the uj (dot-dashed line) 
versus the nucleon momentum for kp = 225 MeV/c. Note the flatness of the 
total (a + uj) contribution. 

Fig. 5 The free (dotted) and the HF charge responses at q — 300, 500, 800 and 1000 
MeV/c. The HF response is displayed in three versions: with the true self- 
energy (p|-[T3|) (solid), with the parabolic self-energy ([27|,[2!|) (dashed) and with 
the analytic method discussed in the text (dot-dashed). Note how close the HF 
curves are, but for the q = 300 MeV/c case at very low energy. The calculation 
is performed for kp = 225 MeV/c, which should roughly correspond to 12 C. 

Fig. 6 The sum rule S L , the mean energy X L and the variance a L of the charge 
response in HF (crosses), TD (squares) and TD on a HF basis (diamonds) 
compared to the free RFG results (dashed line). Note the expected fulfillment 
of S L in all instances, the hardening of the charge response induced by the HF 
field but not by the TD correlations (although, in principle, these do not have 
to fulfill the energy-weighted sum rule) and the near coincidence of the variance 
of the RFG with and without TD correlations. In contrast, the HF variance is 
significantly larger than the RFG one, but for the lowest momenta, although 
the difference between the two is slightly reduced by the TD correlations. 

Fig. 7 The two Feynman diagrams contributing in first order to the RPA series: (a) 
ring (direct), (b) exchange. 

Fig. 8 The momentum behaviour of the direct ph interaction of the Bonn potential 
(the small p meson contribution is neglected). 

Fig. 9 The charge response in the HF and TD approximations at q — 300 MeV/c 
and kp = 225 MeV/c. The dotted lines represent the free RFG response. 
Moreover, we display: in (a) Rf F (solid); in (b) R^ ee (dash) and R™ g ' HF 
(solid); in (c) R e ^ ch ' ree (dash) and ^ xch > HF (solid); in (d) i?™' free (dash) and 
i?™' HF (solid). 
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Fig. 10 The charge response in the HF and TD approximations at q — 500 MeV/c 
and kp = 225 MeV/c. The meaning of the lines is explained in Fig. 

Fig. 11 The charge response in the HF and TD approximations at q — 800 MeV/c 
and kp = 225 MeV/c. The meaning of the lines is explained in Fig. ^. 

Fig. 12 The charge response in the HF and TD approximations at q — 1 GeV/c and 
kp = 225 MeV/c. The meaning of the lines is explained in Fig. || 

Fig. 13 The charge response in the TD approximation at q = 300 MeV/c and kp = 225 
MeV/c. In (a): the free RFG response (dotted) and ^ xch ' free w ith the full 
interaction (solid), with o + oj (dashed), with o + uj + p (dot-dashed). In (b): 
.RTD,free with the full interaction (solid) and with the full interaction minus the 
contribution arising from the interference of the pion with the other mesons 
(dashed); also shown is the contribution of the pion alone (dotted), the global 
contribution of a, u and p (long-dashed) and the interference of the former 
with the latter (dot-dashed). 

Fig. 14 ^ xch ' free (solid) and the free RFG response (dotted) in the isoscalar (a) and 
isovector (b) channels at q — 1 GeV/c and kp = 225 MeV/c. Also shown 
are the exchange responses with only the ir (dashed) and with o + u + p 
(dot-dashed). 

Fig. 15 The charge response in the HF and RPA approximations at q — 300 MeV/c 
and kp = 225 MeV/c. The dotted lines represent the free RFG response. 
Moreover, we display: in (a) i?^ F (solid); in (b) /^ ing ' free (dash) and i?£ ing ' HF 
(solid); in (c) ^ xch ' free (dash) and ^ xch,HF (solid); in (d) ^ PA ' free (dash) and 
i?f A ' HF (solid). 

Fig. 16 The charge response in the HF and RPA approximations at q — 500 MeV/c 
and kp = 225 MeV/c. The meaning of the lines is explained in Fig. 

Fig. 17 The charge response in the HF and RPA approximations at q — 800 MeV/c 
and kp = 225 MeV/c. The meaning of the lines is explained in Fig. 

Fig. 18 The charge response in the HF and RPA approximations at q = 1 GeV/c and 
kp = 225 MeV/c. The meaning of the lines is explained in Fig. || 

Fig. 19 The solution of eq. (|55|) at q — 300 MeV/c; the curves correspond to Ren^ D 
(dashed), Ren° sc (dotted) and Ren^ PA (solid). 

Fig. 20 The sum rule S L , the mean energy X L and the variance a L of the charge 
response in HF (crosses), RPA (squares) and RPA on a HF basis (diamonds) 
compared to the free RFG results (dashed line). 
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Fig. 21 The charge response in the RPA approximations at q — 300 MeV/c and kp = 
225 MeV/c. For the meaning of the lines see panel (b) of Fig. |6|. 

Fig. 22 The sum rule S L and the mean energy X L of the charge response in HF- 
RPA (diamonds) and HF-TD (crosses) for k F = 196, 201, 206, 208 MeV/c at 
q = 300, 500, 800, 1000 MeV/c respectively. 

Fig. 23 The shift of the quasielastic peak of the reduced charge response in HF-RPA 
(solid line) and HF-TD (dashed line). The thin curves correspond to kp=225 
MeV/c, the thick curves to the variable kp. 

Fig. 24 The free (dotted) and the HF-RPA (dashed) charge responses corresponding 
to kp = 225 MeV/c are displayed as functions of the energy transfer. The 
solid curves represent the HF-RPA response functions where the values of kp 
are 196 , 201, 206 and 208 MeV/c for q = 300, 500, 800 and 1000 MeV/c 
respectively. 

Fig. 25 ^ PA ' free a t g = 300 MeV/c with the full interaction (solid), with only the n 
(dashed) and with only the a (dot-dashed). In panel (a) kp = 225 MeV/c and 
in panel (b) kp = 196 MeV/c. 
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